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Dynamical Coherent-Potential Approximation (CPA) to correlated electrons has been ex- 
tended to a system with realistic Hamiltonian which consists of the first-principles tight-binding 
Linear Mufhntin Orbital (LMTO) bands and intraatomic Coulomb interactions. Thermody- 
namic potential and self-consistent equations for Green function are obtained on the basis of 
the functional integral method and the harmonic approximation which neglects the mode-mode 
couplings between the dynamical potentials with different frequency. Numerical calculations 
have been performed for Fe and Ni within the 2nd-order dynamical corrections to the static 
approximation. The band narrowing of the quasiparticle states and the 6 eV satellite are ob- 
tained for Ni at finite temperatures. The theory leads to the Curie- Weiss law for both Fe and 
Ni. Calculated effective Bohr magneton numbers are 3.0 ^b for Fe and 1.2 /^b for Ni, explain- 
ing the experimental data. But calculated Curie temperatures are 2020 K for Fe and 1260 K 
for Ni, being still overestimated by a factor of two as compared with the experimental ones. 
Dynamical effects on electronic and magnetic properties are discussed by comparing with those 
in the static approximations. 

KEYWORDS: dynamical CPA, correlated electrons, tight-binding Linear Muffintin Orbital method, iron, 
nickel. Curie temperature, effective Bohr magneton number, quasiparticle states 



1. Introduction 

Understanding of electronic and magnetic properties 
of the system with intermediate strength of Coulomb in- 
teractions has been a challenging problem over half a 
century in condensed matter physics, because simple the- 
oretical approaches are not applicable to the system in 
spite of the fact that many intriguing phenomena are 
found there. ^^^ Iron and nickel are considered to be an 
example of such systems. These metals in fact show the 
properties of both the weakly- and the strongly- corre- 
lated electrons. Photoemission data, for example, show 
the existence of metallic d bands^'^ and the Sommer- 
feld coefficients in the T-linear specific heats show rather 
large values (5-7 mJ/K^mol) as compared with those of 
the noble metal systems.^ The quasiparticle band widths 
however are found to be narrower than the results of 
usual band calculations and a satellite peak is observed 
at 6 eV below the Fermi level in Ni,^ which are not able to 
be explained by a simple band theory,^ suggesting rather 
strong electron correlations in these systems. The same 
features are found in the magnetic properties. Noninte- 
ger values of the ground-state magnetization in Fe and 
Ni are well explained by the band theory, ^° while the 
paramagnetic susceptibilities of these systems follow the 
Curie- Weiss law and their effective Bohr magneton num- 
bers are close to those expected from the local-moment 
model. ^ Large specific heats near the Curie temperature 
Tq are also well explained by the same model. 

The magnetic and electronic properties in the inter- 
mediate regime of Coulomb interactions have been tra- 
ditionally explained by interpolation theories between 
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the weak and strong Coulomb interaction limits. Cyrot^^ 
proposed an interpolation theory on the basis of the func- 
tional integral method which transforms an interacting 
electron system into an independent electron system with 
time-dependent random charge and exchange fields. He 
showed that the static and saddle-point approximations 
to the functional integral scheme can explain the local- 
moment vs. itinerant behavior in magnetism of transition 
metals, as well as the metal-insulator transition. 

Hubbard^^ and Hasegawa^'^ independently developed 
a single-site spin fluctuation theory. They adopted a high 
temperature approximation (i.e., the static approxima- 
tion), and treated the random charge and exchange po- 
tentials by making use of the coherent potential approx- 
imation (CPA). The theory qualitatively described the 
magnetization vs. temperature curves, the Curie tem- 
perature, as well as the Curie- Weiss susceptibility in Fe 
and Ni. The theory, however, reduces to the Hartree-Fock 
approximation at zero temperature because it relies on 
the static approximation. This means that the theory 
does not take into account the ground-state electron cor- 
relations as discussed by Gutzwiller,^'' Hubbard, ^^ and 
Kanamori.^^ Furthermore, the quasiparticle bands and 
the satellite peak do not appear in the theory using the 
static approximation. 

We proposed the dynamical CPA which fully takes 
into account the electron correlations within the single- 
site approximation, and clarified the qualitative fea- 
tures of dynamical effects using a Monte-Carlo sampling 
method. ^^ More recently, we developed analytic method 
to the dynamical CPA,^^ adopting the harmonic approx- 
imation (HA).^^ The latter is based on the neglect of 
the mode-mode couplings between dynamical potentials 
in solving an impurity problem in an effective medium. 
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The HA interpolates between the weak Coulomb inter- 
action limit and the atomic limit. Especially it describes 
the Kondo behavior quantitatively in the strong corre- 
lation limit. ^° We showed within the single band model 
that the dynamical CPA-I-HA yields the band narrowing 
of quasiparticle states and the satellite peak in Fe and 
Ni, which were not explained by the early theories with 
use of the static approximation. The theory was however 
based on the single-band Hubbard model. Quantitative 
calculations of transition metals and alloys with use of 
the realistic Hamiltonian has not yet been made even 
within the single-site approximation. 

In the present paper, we extend the dynamical CPA 
to the multi-band case adopting the first-principles 
tight-binding (TB) linear muffintin orbitals (LMTO) 
method. ^^'^^ The modern band theory is based on the 
density functional theory (DFT) which allows us to ex- 
press the ground-state energy as a functional of the spin 
and charge densities of the system. ^'^ In the local density 
approximation (LDA),^** one approximates the energy 
functional with the energy function of the density. The 
LDA exchange-correlation potentials obtained from the 
electron gas system have much simplified the electronic 
band-structure calculations in solids. The TB-LMTO 
method allows us to construct the first-principles tight- 
binding one electron Hamiltonian, and to calculate the 
LDA band structure. We adopt the TB-LMTO Hamil- 
tonian to describe the noninteracting part of the Hamil- 
tonian, and take into account the intraatomic Coulomb 
and exchange interactions between d electrons which are 
dominant among electron-electron interactions. 

Similar theoretical approach has been developed in the 
problem of the metal-insulator transition in infinite di- 
mensions.'^''' The approach called the dynamical mean 
field theory (DMFT) is equivalent to the dynamical CPA, 
as we have shown recently. ^^ The present theory there- 
fore should be equivalent in principle to the DMFT com- 
bined with the LDA-I-U scheme in the band theory. ^^ 
The merits of the present approach may be summarized 
as follows. (1) The dynamical CPA can treat the trans- 
verse spin fluctuations for arbitrary d electron number 
at finite temperatures , while the standard DMFT com- 
bined with the quantum Monte-Carlo method (QMC) 
cannot treat them because it is based on the Ising-type 
Hubbard-Stratonovich transformation.^® Because of the 
reason, the DMFT calculations for Fe and Ni have been 
performed so far without taking into account the trans- 
verse spin fiuctuations at finite temperatures.^^ (2) The 
HA which we adopted to solve the impurity problem is 
an analytic approach from the high-temperature limit. 
The approach is suitable for understanding the finite- 
temperature magnetism because the zero-th approxima- 
tion to the HA describes the magnetic properties much 
better than the Hartree-Fock one. There is no corre- 
sponding approach in the DMFT. (3) Because of the an- 
alytic theory, we can calculate the excitation spectra up 
to the temperatures much lower than those calculated by 
the QMC, using the Padc numerical analytic continua- 
tion method. "^^ 

In the following section, we introduce a TB-LMTO 
Hamiltonian with intraatomic Coulomb interactions. In 



§3, we formulate the dynamical CPA to the realistic 
Hamiltonian on the basis of the functional integral tech- 
nique.'^^ Applying a generalized Hubbard-Stratonovich 
transformation"^^ to the free energy, we transform the 
interacting electrons into an independent electron sys- 
tem with time dependent random fields. Introducing an 
effective medium into the time dependent Hamiltonian, 
we will make a single-site approximation. We determine 
the medium solving a self-consistent equation, called the 
CPA equation. "^^ In §4, we adopt the HA to calculate 
the dynamical part of the free energy, and derive the 
analytic expressions of the free energy, the dynamical 
CPA equation, and other thermodynamic quantities. In 
§5, we present the numerical results of calculation for 
Fe and Ni. The calculations have been performed by us- 
ing the second-order dynamical CPA (i.e., the dynam- 
ical CPA-I-HA within the second-order dynamical cor- 
rections). We explain the band narrowing of the quasi- 
particle states, the incoherent satellite peak at 6 eV 
below the Fermi level in Ni. We also present the re- 
sults of calculations for the magnetization vs. tempera- 
ture curve, the paramagnetic susceptibility following the 
Curie- Weiss law, and the amplitude of local moments. 
We clarify the quantitative aspects of the theory com- 
paring with the experimental data, and examine the dy- 
namical effects on various quantities comparing the dy- 
namical results with those in the static approximation. 
The last section 6 is devoted to summarize the dynami- 
cal CPA and TB-LMTO Hamiltonian approach, as well 
as the dynamical effects in Fe and Ni. 



2. TB-LMTO Hamiltonian 

We adopt in the present paper the first-principles TB- 
LMTO method^^ to construct a realistic many-body 
Hamiltonian. In this case, atomic basis function with 
orbital L on site i, XiLii" ^ Ri), are constructed from 
a muffintin atomic orbital ipiLir ^ Ri) on site i with 
an atomic level E^^iL, and a tail function outside the 
muffintin potential h"j^,^j^ as 

X^L{r-R^ = ^^L{r-R^) + J2^7L'{r-RJ)h^L'^L, (1) 



"PtLir - Rt) = ^iL{r - Ri) + ^Pihir - Ri)o'^L 



(2) 



Here the wave function in the interstitial region has 
been neglected because of the atomic sphere approx- 
imation. The atomic wave function ipiLir — Ri) and 
its energy derivative ipiL{r) are defined by ipiL{r) — 
(l)iL{E,yiL,r)YL{r) and <ftL{r) = <j)iLiE^iL,r)YL{r), 
where Yl^t) is the cubic harmonics with L = (l,m), I 
being the azimuthal quantum number and m being an 
orbital index for I. (f>iL(E,r) is obtained by solving the 
radial Schrodinger equation with energy E. The energy 
Ei,iL is chosen to be the center of gravity below the Fermi 
level for each orbital. The atomic orbitals {ipiL} are nor- 
malized in the atomic sphere as {(piL\^iL') = Sll'- The 
tail coefficients h%,^j^ in eq. (1) are determined in such a 
way that the orbital XiL is continuous and differentiable 
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on the sphere boundary at each sphere. The cocfhcient 
o"^ in eq. (2) is determined so that the orbital XiL is well 
localized. We adopt here the nearly orthogonal represen- 
tation (i.e. ofj^ =0), in which the orbitals Xjl's become 
orthogonal up to second order in hiLju ■ The TB-LMTO 
Hamiltonian matrix is then written as 

H.LjV = (x»L|(-V2+«(r))|x,L> 

= ^iLSij^LL' +tiLjL>. (3) 

Here v{r) is a LDA potential, e^^ is an atomic level, and 
tiLjL' is a transfer integral between orbitals XiL and XjL' ■ 
When we adopt the density functional theory, the one- 
electron Hamiltonian (3), especially the atomic level e^L, 
contains the effects of strong intratomic Coulomb inter- 
actions in general. According to the LDAH-U interpreta- 
tion by Anisimov et. al.,^^ the atomic level e"^ for non- 
interacting system is obtained from the relation. 



dE] 



LDA 



dE^^BA 



HL 



OUiLa OUiLa 

Here ruLcr is the charge density at the ground state, 
-Elda is the ground-state energy in the LDA, and -Elda 
is a LDA functional to the intraatomic Coulomb inter- 
actions. Among various forms of £^ldA' ^^ adopt the 
Hartree-Fock type form"^^ since we consider here an itin- 
erant electron system where the ratio of the Coulomb 
interaction to the d band width is not larger than one. 



E, 



LDA 



j mm' a 

+ 9X111 ^(U-^hjdrijd 



(5) 



Here rijd = X]m<T'^i'™'^/2(2' + l) with I = 2. U and J are 
the orbital-averaged Coulomb and exchange interactions 
defined by 

1 



U 



{U-J) = 



{21 + 1)2 ^, 



U„ 



biS''^" 



21(21 + 1) 



Ju 



(6) 
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where Umm' and Jmm' are orbital dependent intraatomic 
Coulomb and exchange integrals for d electrons. From 
eqs. (4) and (5), we obtain the atomic level e°^ for non- 
interacting system as 



ejL 
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ndS, 



12 ■ 



(8) 



(2/ + 1), 

Note that Ud denotes the total d electron number per 
atom. 

The Hamiltonian which we consider here can be writ- 
ten as 

H = Ho + Hi. (9) 

The tight-binding Hamiltonian for noninteracting system 



Ho is given by 

^0 = J2^^"^L - A*) n^La + Yl ^^^1^' "'Iha^jL'a ■ (10) 
iLa iLjL'a 

Here we have introduced the chemical potential ^ for the 
calculation of the free energy. a\^^ [aiLa] is the creation 
(annihilation) operator for an electron with orbital L and 



spin a on site i, and fiiLa 



^iLa 



UiLa- is a charge density 



operator for electrons with orbital L and spin a on site 
i. We have neglected the change of the transfer integrals 
due to electron-electron interactions. 

The interacting part Hi in eq. (9) consists of the in- 
traatomic Coulomb interactions between d electrons. 



-ffi = ^ I y^L/pn^imT 



'^ilml 



1 



^ 2^ (Ui - -:J)niimnii„i' 



E 

m'>mi' 



Js. 



ilm ' ^ilm' 



(11) 



Here Uq (Ui) and J are the intra-orbital (inter-orbital) 
Coulomb interaction and the exchange interaction, re- 
spectively, fiiim (siim) with Z = 2 is the charge (spin) 
density operator for d electrons on site i and orbital 
m, which is defined by fium = J2a '^Uma {Siim = 
Sa-Y ALai'^l'^)a-iO,tL-t), CT bciug thc Pauh spin matrices. 



3. Functional Integral Approach and Dynamical 
CPA 

Thermodynamic properties of the system are calcu- 
lated from the partition function, which is given by 



Z = Tr 



Tcxp 



/3 



{Ho{t) + Hi{t)) 



(12) 



Here (3 is the inverse temperature, T denotes the 
time-ordered product (T-product) for operators. Hq{t) 
{Hi{t)) is the interaction representation of Hamiltonian 
Ho {Hi). 

The functional integral method is based on a Gaussian 
formula for the Bose-type operators {ft^}. 

/ ^ (^m-^mm' O^m' 
g mm' ^ 



detA 



rM 



m 
/^\'^m-^mmi''^mi' ^^ ^G^mi-^mim' 



XjYi' j 



.(13) 

Here Amm' is a M x M matrix, and {a;™} are auxiliary 
field variables. Discretizing the integral with respect to 
time in eq. (12), and applying thc formula (13) to the 
bosc-type operators at each time under thc T-product, 
we obtain a functional integral form of thc free energy T 
as 



-fir 



■ N 21 + 1 

n n ^i^m{r)5Crm{T) 



_z^l m— 1 



Z\i{r)X{T)) 
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i mm' '^ a \ m' / 



W'.(24) 



xyz 



/ , <iima(T)-Dj^,ui^il;im'a(T, 



, (14) 



The functional integrals in the Fourier representation in 
eq. (21) is given by 



^°(l(T),C(T))=Tr Texp 



-j H{T,t^K)dr j, 



(15) 



2/+1 



[ n '^c. 



n r-n^ n ^c.™(o) 



(47r) 



m— 1 



■ff (t, ^ -«C) =X! [(^»'i^/^ ^2^ i^imm' C,im' {t)5i2J fliL (? 
iL m' 



11 (4^)2^ + 1 ^ C.m(*C^n) 
.n=l ^ ' 



.(25) 



(16) 



iLjL' a 

Here N is the number of sites, rhiL — 2siL, and 
CimiT) (^i„i (''')) is 8-1^ auxiliary field being conjugate with 
iniL{T) (rhiL{T)) for 1 — 2. The functional integrals in 
eq. (14) are, for example, defined by 



pi+i 



ll^'5Cim(7 



jri— 1 



N' 






(4^) 



/N' 



, (17) 



where 2Z + 1 in the square roots denotes the number of d 
orbitals (i.e, 21 + 1 = 5). detA^ is the determinant of the 
(2/ + 1) X {21 + 1) matrix Aimm' ■ Tn denotes the n-th time 
when the time interval [0, (3] is divided into N' segments. 
The matrices Aimm' and Bf^^, [a ~ x, y, z) are defined 
as 



Here the field variable C,im{i^n) {£,ima{i^n)) denotes 
the n- frequency component of Cim(T) {^ima{T)), and 
d'^Cimiii^n) = dRe(im{iuJn)dlm(imiiLu„). The energy 
functional E[^, C] in eq. (22) consists of the noninter- 
acting term (the first term at the r.h.s. (the right-hand- 
side)), the scattering term due to dynamical potential 
(the second term), and the Gaussian term (the third 
term). Sp in the second term at the r.h.s. of eq. (22) 
means a trace over site, orbital, frequency, and spin, g in 
the second term denotes the temperature Green function 
for noninteracting system Hq- The dynamical potential 
V is defined by eqs. (23) and (24), and <Ta in cq. (24) 
denotes the a component of the Pauli spin matrices. 

In the effective medium approach, ^^ we introduce a 
coherent potential 



i^)iLncrjL'n'cr' — ^ Lcr{i^n)SijS LL' Snu'^a 



(26) 



into the energy functional -E[^,C]i and expand it with 
respect to w — S as 



E[tC]=^ + yEd^„Q + ^E 






(27) 






Bl 



UaSmm' + '^(l ~ 0mm') 



(18) 
(19) 
(20) 



Equation (15) is a partition function for a time- 
dependent Hamiltonian H{t, ^, —iQ of an independent 
particle system. Note that we have introduced a mag- 
netic field hfm for convenience. 

In the Matsubara frequency representation, the free 
energy !F is written as 

» N 21+1 

e"^^ = / [ n n ^^jmSQm] exp [-mt C] ] , (21) 
•^ j = lm=l 



-\T.T. 



^1771 y.'^^'n ) -^imm' <,ni' \'^^n) 



(IUJ„ 



(22) 



+ 2_^ ?ima(*^«)-°imm'Mm'Q(* 
a 

iv)2Ln<jjL'n'a' ^ VjLcra'ii^n- i^n')Sij6LL' , (23) 

ViLaa'{iuJn) = -x2j*^»™ni'Cim'(»W„)(5;2<5<T<T' 



Here the zero-th order term JF is a coherent part of the 
free energy which is defined by 

T = -l3-'^\nTr{c-f^"°) - /3-iSpln(l - Sg) . (28) 

Note that the coherent part does not depend on the dy- 
namical potential. 

The next term in eq. (27) consists of a sum of the 
single-site energies Ei [^^ , C,.i\ , which are defined by 

+ 4 E E [Qm{i^n)A^mm'C,im'{i'^n) 
n mm' 
xyz 

+ Yi*ma{i^n)BZnm'£.im'a{iUJn)] ■ (29) 
a 

Here tr means a trace over orbital, frequency, and spin on 
site i. Svi = Vi — Sj, and Vi (Si) is the dynamical (coher- 
ent) potential on site i. F, is the site-diagonal component 
of the coherent Green function defined by 

{Fi)jLnaj'L'n'a' ~ FiLa{i^n)SijSij' 5lL' 5nn' Sera' , (30) 
FiLa{i^n) = [{9 — S) ]iLnaiLna ■ (31) 

The last term in eq. (27) denotes the higher order 
terms in expansion. 



AE^~(3'^Sp\n{l-tF') 



(32) 
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(33) 4. Harmonic Approximation to the Dynamical 
CPA 



and F' is the off-diagonal coherent Green Function de- 
fined by We can rewrite the free energy (35) by means of an ef 

fective potential projected onto the zero frequency vari- 



{F )iLnajL'n'cr'— [{g — S) JiLncrjL'n'cr' (1— % )(5crcr' • (34) 

The dynamical CPA is a single-site approximation 
which neglects the intersite dynamical correlations AE. 
The free energy is then written as 

^CPA =.^ -E /5"' In l\ll^^rmSQm]e-(^^^l^^^'^^^ . (35) 

The dynamical coherent potential ^^^^(JtiJn) should be 
determined so that the nonlocal corrections AE vanish 
in average. This means that 



ables ^,„ == €jm(0) and Q^ = Gm(0). 

■^CPA = ^ 



11 V (A^Vl+1 ii'^^" 



(4^) 



f3^^+^dctA 



;2/+l 



Y[dC„ 



-PEitC) 



(44) 



where 



((-)> 



iu) = , 



Note that we have redefined Tcpa and T by those per 
(36) site, assuming that all the sites are equivalent to each 
other. Furthermore we omit here and in the following all 
the site indices for simplicity. 

The effective potential E {$,,() in eq. (44) consists 
of the static part Est{^,C) and the dynamical one 

-£^dyn(l, 0- 

(45) 



(37) 



The above condition called the CPA equation is writ- 
ten as 



('^iLa(*^")) = FiLa{i^n) 



.(«) 



G^L(^^«) - [iFf 



Svi 



liLnaiLna 



(38) 
(39) 



Here the l.h.s. (left-hand-side) of eq. (38) is a tempera- 
ture Green function for an impurity system in the effec- 
tive medium, whose Hamiltonian is given as follows. 



i^), 



H^^'ir) =. H{t)+H'^>{t) 



- f dr'J2 4LA^)^^Ur - T')a,UT'), (40) 

•^0 La 

iL<7 



iLjL' <7 

rl3 



jLcr 



E/ dT'a]^jT)^,Ur-r')a,Ur'), (41) 



Hi (t) = y^ Up niim-\ {T)hil.mi (t) 
m 

'm'>m' 



(42) 



It should be noted that the CPA equation (38) is equiv- 
alent to the following stationary condition. 



(5I]iLcr(«W„) 
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(43) 



Est{i,0 = -/?"'trln[l-r5«oi^,] 

xyz 



Q-PEdyn (^, = 



11 (27r)2'+i ^^'""('^") 
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J](i^Cm(«t^«) 



/3 



E E (Cm(*^n)^mm'Cm'(«^n) 
n^^O rrnn' 



xD exp 

4 ■' — ' ■' — ' \ 

•- n^O mm' 

D = detlSnn'SLL'Saa' 

-~^VLaa"{iuJn-ii^n')9La"L'a'{i^n')j ■ (48) 



Here dvp in eq. (46) is defined by 5vo — u(0) — S 

{5va)LnaL'n'a' = (^Lcrcr' (0) - Y.La{iLOn)5aa')6LL'5n 



(49) 



VLa-a' (0) is thc static potential, while v in eq. (48) is the 
dynamical potential without zero frequency part. 

VLaa'{it^n - ii^n') = VLaa'{ii^n - ii^n') - VLaa'{0)Snn' ■ (50) 

Furthermore, cjLaL'a'ii^n) in cq. (48) is the Green func- 
tion in the static approximation defined by 



iLnaL'na' 



(51) 



where the coherent Green function F is defined by eqs. 
(30) and (31). 

In the functional integral approach, we first have to 
calculate the determinant (48), and second have to eval- 
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uate the functional integral in eq. (47). In order to im- 
plement these calculations, we expand the determinant 
(48) with respect to the frequency modes of the dynam- 
ical potential fLo-cr' as follows. 



'XgLa"a"'{n- v)vLo 
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D = l + Y,{D.-l) 

+ Y^ (A.,.' -D,- Ay + 1) 



(52) 



D^ = det 



^LL'&aa'^nn' - '^^ (vhaa" (j'^v)^n- 



+VLaa"ii^-i^)Sn-n'-u)gLa"L'a'ii^n') , (53) 



We assumed in the above expression that the orbitals {L} 
form an irreducible representation of the point group of 
the system, so that gLaL'a'{ii^n) = gLaa'iiuJn)SLL' (see 
eq. (51)). Furthermore here and in the following, we write 
the frequency dependence, for example, of gL^^iliujn) as 
gLaa' (n) for simplicity. 

The determinant D^{k,m) is expanded with respect 
to the dynamical potentials as follows. 

D^{k,m) = l + Di^\k,m) + D^^\k,m) + --- , (59) 



Q:i7i ■■■OLnIn 
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det 
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^aa' ^nn' 



-^{vLaa" {it^i)5n-^'.y+VLaa" {it^-i^5n-^i'-v)gLa"L'cr' («^n') 
a" 



X Va„ [v, m)u^„ {~v, m)D^^J^y {v, k, m) . (60) 

Here the subscripts a^ and 7^ take 4 values 0, x, y, and 
z, and 



Wo 



(j/, m) = --iy^^Amm'Cm'{l^)Sl2 



(61) 



+VLaa"{i0J-,.')5n-^'.-,.')gLa"L'a'{iuJn')]- (54) ""('^''™) 



-7;YB^nm'S.rn'aMSl2 , {a^X,y,z). (62) 



The first term at the r.h.s. of eq. (52) corresponds to 
the zero-th approximation (i.e. the static approximation) 
which neglects dynamical potentials. The second term 
is a superposition of the independent scattering terms 
of dynamical potential VL^^i{iLU^). Higher order terms 
describe dynamical mode-mode couplings. 

We adopt here the harmonic approximation^^ which 
neglects the mode- mode coupling terms in eq. (52). We 
have then 



Note that the subscript {aj} in eq. (60) denotes a set 



of (ai7i • • • a„7„). The expressions of D/^ , (:/, fc, to) are 



given in Appendix A. 

Substituting eq. (59) into eq. (56) and taking the Gaus- 
sian average (47), we have 



21 + 1 



D^ 



i^dy„(^,C)--r'ln 



l + Y,{D,-l) 



(55) 



51 51 ^ [ nn^"'(''''^)^'>'»(^^'' 

21+1 u-1 



X 



The approximation yields the result of the second-order 
perturbation in the weak Coulomb interaction limit, and 
describes the Kondo anomaly in the strong interaction 
limit. 2" 

Let us now calculate £)^ in eq. (55). The determinant 
Di, in the harmonic approximation is written by a prod- 
uct of those of the tridiagonal-type matrices as 

v-l [21 + 1 



[nn^{Sr"('^'^'"^)]-(63) 

m—l k—0 

Here {l{k, ■m)}{k = 0, • • • ,iy — l,m — 1, • • • , 2Z -I- 1) are 
zero or positive integer, satisfying X]fcm'(^''^) ~ "• 
Y[^Va.{h',m)v^.{—i>,'m) are the products of 
Vai{i',m)v^.{—i>,m) belonging to the m-th orbital 
block. Calculations of the Gaussian average of the 
dynamical potentials are given in Appendix B, and we 
reach the following expression. 



D,^ 



n 



l[D,{k,m) 



k=0 im=l 



(56) 



Diyikjiri)- 



1 1 

a-,^(y,m) 1 1 

afc(i^,TO) 1 1 

a,^fc (i/,to) 1 1 

a2^+fc(j^,m) 



D 



in) 



.(57) 



A. = 


1 


(2/3)", 


21+1 iJ-1 ' 

nil 

m=l k=0 . 



l + D 



{1) 



D 



■(2) 



(64) 



E E E 

Efcm l{k.m)^n {oij (fc,m)} P 



l{k,'. 



n ^™€'™^p£t'"i\(-'^'-) 



i=i 



.(65) 



Here 1 in the determinant is the 2x2 unit matrix, 
a„ (y, m) is a 2 x 2 matrix defined by 



Here j denotes the j-th member of the (fc, m) 
block. P denotes a permutation of a set {(j, fc,TO,)}: 
P{(j, fc,TO-)} = {(jp, fcp,TO-p)}, a-p-i means an rearrange- 
ment of {aj{k, to)} according to the inverse permutation 
P^^. Note that aj{k,m) takes 4 values 0, x, y, and z. 
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(66) 



Equations (55) and (64) determine the dynamical poten- 
tial ^dyn(€, 0- 

The free energy (44) is written alternatively as 



^CPA -^-/S^^ln 



^ /j^'+ideti?" -rr ^^ 



(4^) 



-/3i?cff(^) 



(67) 



In the itinerant electron system, spin fluctuations plays 
an important role, and we may neglect the thermal 
charge fluctuations making use of the saddle-point ap- 
proximation to the static charge flelds Cm • We have then 
Ecsii) = E{^,C*)- The saddle point value Cm is deter- 
mined from dE{^, C*)/dCm = 0: 



■<m =^i(0 =X1^^'^(^) 



■n-L, 



■iO^^T^GUr 



(68) 



(69) 



In order to reduce the number of variables, we ne- 
glect the out-of-phase thermal spin fluctuations between 
different orbitals on a site, and take into account their 
in-phase fluctuations. This can be made by introduc- 
ing a large variable S.a = J^m^rna- Inserting 1 = 

I [Ha d£.adXa] exp[-27riAa (^a - J^m ^rna)] mtO Cq. (67), 

and replacing variables ^ma with ^a/{'2l + 1) in the non- 
Gaussian terms of Ecs{^), we reach 



^CPA = ^-rMn| [n^^de.]e-^^<=«'(^), (70) 

EeS (€) = E,t (^) + i?dy„ (^) , (71) 

^^t(€) = -4 E ^n[{l-5v^{0)F^{n)){l-6vLi{0)FLi{n)) 



ljie^FL^{n)FLi{n) 



-{Uo-2Ui + J)J2nLii) 



-(2C/i - J)M^r + JUi + JtAt ■ (72) 



Here Jx 

Uo/{2i ^ 

VLaiO) : 



J-u 



'y^- J± = (1 - 1/(2/ + 1))J, J. - 

1) + J_i, SvLa{0) = wl<t(0) - ^Lcrin), and 
vo{0,m) + avz{0,m). The charge densities, 
%(!) and ni{$,) are deflned by %(!) = J^a^Lai^) and 
^i(^) = Z]m"i(^)- Furthermore -Edyn (I) is given by 
eq. (55) in which Cm (Cma) has been replaced by ifiLi^) 

(Ca/(2/+l)). 

The CPA equation in the HA is obtained from the 
stationary condition (43) with the free energy (70). 



(GLain)) = FLain) 



(73) 



(GLa{n)) = (gLaain) TTI^ f \ 



dyn 



(74) 



Here KLa(n) = 1 - FLa{ny^HLa{n) and HLa{n) = 
5FLa{n) / 5YiLa{n). The average (^) at the r.h.s. of eq. 
(74) is now deflned by a classical average with respect to 
the effective potential (71). 



[n^^"](^ 



-PE,s{£.) 



[\{d^. 



-PE^siO 



(75) 



Substituting eq. (55) into Eq. (74), we obtain the expres- 
sion 



I ? 

{GLcr{n)) = ( gLaa{n) 



KLa{n)ST.La{n) 



i + E.(^.-i) 



(76) 



The local charge and magnetic moment are obtained 
from d!FcPA/de1 and —dTcvk/dhf^. Making use of the 
stationary conditions of J^cpa with respect to Cm ^^d 
Sicr, and the CPA equation (73), we reach 



{^L 



^E^-(^ 



/3 



<^i) = ^E'^^^-(^ 



(77) 



(78) 



In particular, the I = 2 components of local charge and 
magnetic moment are expressed as 



{ni) = {nim , 



(79) 



(rhi) = (0 . (80) 

The amplitude of charge and local moments for d elec- 
trons are calculated from the formulae. 



(»f) = (n,)+2i:^^^'* 



Y. 



311 ^-^ 311 

mm' 



(81) 



3(n,)-6E 



(y^mmi 



, \^^mrn' ^'-'mni' / 



(82) 



Here we have introduced for convenience orbital- 
dependent Coulomb and exchange interactions Umm' and 
Jmm' into the interaction Hi to derive the expressions. 
Making use of the stationary conditions of J^cpa and in- 
tegrations by parts, we obtain 



(n?) 



("-; 



(0) + 1 E("i(0'> + E'("i(o%'(€)> 



2{2l + l)y'''' i3J^ 



2E 



dE, 



dyn 



dU„ 



E' 



dE, 



dyn 



dU„ 



(83) 
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(rhf) = -Sifni^)) - -J2{hUe 



2(21 + 1) y^'^ I3J, 



-6E 



2/ 

dEdyn 



a^x.y 



E 



dEdyn 



+2 



811 

_ ^'-^ mm 

dEdyn 
yy u rn.m ' 



(84) 



Here [dE^yn/dUmm'] means taking derivative of -Edyn 
with respect to Umm' fixing the static potentials wlo-ct' (0). 
In the HA, these values are obtained from eq. (55) as 



dE, 



dyn 



dUrr 






'-'^ rn.r. 



(85) 



l + Y,(D.-l) 



BE, 



dyn 



dJ„ 






dJr, 






(86) 



The entropy is calculated from /3^9JFcpa/9/3 as 



'dp 



dp 



^^j[I[^l%d^^--^^EM^^-^''^^^^^^^-lm 



Here 



P- 



P 



'dp N 
dE,ff 



^Spln(g-i-S) + ^^Fi.(r 



(88) 



n La 



dp 



(trln(l-,5«oF))-/3(^dyn) 

"9(/3Sdyn 



p 



dp 



(89) 



The first term at the r.h.s. of eq. (87) {i.e. eq. (88)) 
is the contribution from the coherent free energy and 
reduces to the entropy Sq for noninteracting electrons 
when S^cr — > 0: 

S'o = -2 f dwp^ {uj)[f{uj) In f{uj) 

+ (l-/(c.))ln(l-/(c.))] , (90) 

where p'^{lo) is the total density of states per spin for 
noninteracting electrons, f(uj) is the Fermi distribution 
function. The second term in eq. (87) {i.e. eq. (89)) is 
the entropy due to the temperature dependence of the 
effective potential. [d{PEdyn)/dP]^^ in eq. (89) means to 
take the derivative with respect to /3 fixing the frequency 



iuJn and the coherent potential Sio^(«cij„). It is given in 
the HA as 



T^(") 



P 



5(/3Sdyn) 



dp 



oo 



(91) 



The third and fourth terms in eq. (87) produce the mag- 
netic entropy due to thermal spin fluctuations. 

The thermodynamic energy is obtained from the rela- 
tion {H - fiN) = J^cPA +P~'^S as 



{H - MiV) 



;^EE*^"^^'^'^") 



P 



n Lo 



{Uo - 2Ui + J)J2{fiL{if 



-{2u,- JHM^r) -J2jj{e^) - -^ 



PJa 



CO oo 



EE"^^' 



(«) 



+r ^ { ^^^^ ) .(92) 

\n-^(U.-i)/ 

The first term at the r.h.s. of eq. (92) is the coherent 
contribution of the kinetic energy, the second term corre- 
sponds to the double counting correction in the Hartree- 
Fock energy. The last one is the dynamical correction to 
the energy. 

The sum rule no = ^i^{nL) determines the chemical 
potential for a given valence electron number ng. The 
CPA equation (73) and effective potential (71) with eqs. 
(72), (55), and (64) form the self-consistent equations to 
determine the dynamical coherent potential {T,L^{iLu„)}. 
After having solved the self-consistent equations, we can 
calculate the magnetic moments and charge from eqs. 
(77), (78), and (80), the square of local charge and spin 
fluctuations from eqs. (83) and (84), as well as the other 
thermodynamic quantities (see eqs. (70), (87), and (92)). 



5. Numerical calculations: Fe and Ni 

The simplest approximation to the dynamical CPA is 
to neglect the dynamical potential -Bdyn (C) in the self- 
consistent equations. This is called the static approxi- 
mation and may be justified in the high temperature 
limit. The next approximation is to add the dynami- 
cal potential -Edyn (C) by taking into account the higher- 

order terms D (n > 1) in a series expansion (64). We 
have taken into account the terms up to the second order 
(n < 2) in eq. (64). We call this level of approximation 
the second-order dynamical CPA. Within the approxi- 
mation, we have performed numerical calculations for Fe 
and Ni in order to examine the quantitative aspects of 
the theory and the dynamical effects on their electronic 
and magnetic properties. 
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Fig. 1. Densities of states (DOS) calculated by the LDA and TB- 
LMTO method. Dashed curve: local DOS for Cg electrons, dotted 
curve: local DOS for t2g electrons, solid curve: total DOS con- 
sisting of 4s, 4p, and 3d orbitals 
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Fig. 2. Single particle excitation spectra (DOS) for d electrons in 
the paramagnetic Fe. Results for the LDA, the static approxima- 
tion, and the 2nd-order dynamical CPA are shown by the dotted 
curve, the dashed curve, and the solid curve, respectively. 



We obtained the intraorbital Coulomb interaction Uq, 
interorbital Coulomb interaction Ui , and exchange inter- 
action energy parameter J from the parameters U and 
J in the LDA + U via the relations: Uq = U + 8J/5, 
Ui ^U - 2J/5 and J ^1. (Note that Uq ^ Ui + 2 J.) 
We adopted in the present calculations the LDA+U 
values used by Anisimov et. alP''^ ; U ~ 0.1691 Ry 
and J = 0.0662 Ry for Fe, and U = 0.2205 Ry and 
7 = 0.0662 Ry for Ni. These values yield t/o = 0.2749 Ry, 
C/i = 0.1426 Ry, J = 0.0662 Ry for Fe, and Uo = 0.3263 
Ry, Ui = 0.1940 Ry, and J ^ 0.0662 Ry for Ni, respec- 
tively. 

In the numerical calculations, we adopted an approxi- 
mate expression of the coherent Green function'^® 



FlAu) = 



PL{.f)de 



IUJ„ 



e - 'ELa{iL^n) 



(93) 



The expression takes into account the effect of hybridiza- 
tion between different I blocks in the nonmagnetic state 
via the local densities of states PL(e), but neglects that 
in the spin polarized state. Moreover, we adopted a de- 
coupling approximation^'' to the thermal average of the 
impurity Green function in the dynamical CPA equation 
(73). 



(GL.(n,6,Ci)) = 




,(ei)).(94) 



The approximation is correct up to the second moment 
and reasonably describe the thermal spin fluctuations. 
We have solved the dynamical CPA equation for the bcc 
Fe using the expressions (93) and (94). The densities of 
states (DOS) for 3d, 4s, and Ap states were calculated by 
using von Barth-Hedin LDA potential. The total DOS 
and the d DOS for Cg and t2g electrons are shown in 
Fig. 1. Single-particle excitation spectra have been calcu- 
lated by using the Pade numerical analytic contribution. 
Figure 2 shows the calculated d DOS of paramagnetic 



Fe at T/Tq = 1.19 as the single-particle excitation spec- 
tra. The DOS in the static approximation is broadened as 
compared with the LDA result in the nonmagnetic state 
because of the strong thermal spin fluctuations. The dy- 
namical charge and spin fluctuations produce a satellite 
peak around uj — —0.5 Ry (=—6.8 eV), and suppress the 
band broadening by about 22% as compared with the 
static one. The existence of the satellite peak is consis- 
tent with the previous results of the ground-state calcu- 
lations'^^ as well as those at finite temperatures.^® The 
d band width in the present calculations, though it is 
strongly reduced as compared with the static one, is com- 
parable to that of the LDA calculations, while the XPS 
experiments suggest the 10% reduction of the width as 
compared with the LDA results, and no dip at a; = —0.1 
Ry.'"' These inconsistencies may be attributed to an over- 
estimate of the local exchange splitting above Tc ■ 

Below the Curie temperature, the up and down DOS 
are split as shown in Fig. 3. In the up-spin band, the 
satellite peak at lu — —0.45 Ry remains, and the quasi- 
particle bands at a; « —0.2 Ry shifts to the Fermi level 
as compared with those in the static approximation, 
showing the band narrowing. The satellite peak for the 
down-spin band disappears because of a large value of 
|ImEio^(z)| in this energy region. These behaviors are 
consistent with recent QMC calculations without trans- 
verse spin fluctuations.^® 

It is not easy to calculate the DOS at low tempera- 
tures in the QMC calculations. The present approach al- 
lows us to investigate the DOS even at low temperatures. 
Figure 4 shows the DOS at T/Tc = 0.3. The DOS in the 
static approximation approaches to the Hartrce-Fock one 
with decreasing temperature, but arc still broadened at 
this temperature by thermal spin fluctuations. Dynami- 
cal terms suppress the thermal spin fluctuations and de- 
velops the quasiparticle states, so that sharp peaks of 
Cg electrons appear at oj ^ ±0.15 Ry in the DOS. The 
present calculations reduce to the 2nd-order perturba- 
tion theory at T = 0, so that the DOS in Fig. 4 is close 
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Fig. 3. Up and down d DOS in the ferromagnetic Fe at T/Tq = 
0.7. Results for the static approximation are shown by the dotted 
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Fig. 4. Up and down d DOS in the ferromagnetic Fe at T/Tq 
0.3. 




Fig. 5. Effective potential in the ferromagnetic Fe at T/Tq = 0.5 
on the iix-S,z plane. 




Fig. 6. Dynamical contribution to effective potential in the fer- 
romagnetic Fe at T/Tq = 0.5 on the ^i^-^z plane. 



to those obtained at the zero temperature by Drchal et. 

The effective potential determines the behavior of 
magnetic moments. Figure 5 shows the potential for Fe 
below Tq. It has double minima along z axis, and mono- 

tonically increases with increasing S,±_ = J^"^ + Q . (Note 
that the effective potential is spherical on the xy plane: 
Ecii{S,zi C-l))- The dynamical contribution -Edyn(^) to the 
effective potential is given in Fig. 6. The dynamical part 
shows a 'butterfly' structure; it increases along the z axis 
with increasing the amplitude |^|, while it decreases on 
the xy plane. This implies that the dynamical effects re- 
duce the longitudinal amplitude of spin fluctuations, and 
enhance the transverse spin fluctuations. In fact, we find 
6% reduction of \/{£,1) and 6% enhancement of \/ {S,\) 
at T/Tc = 1.19. 

Magnetic properties of Fe are summarized in Fig. 7. 
Both static and dynamical calculations yield the Curie- 
Weiss susceptibility. Calculated effective Bohr magneton 
numbers are 3.1 ii-q in the static approximation and 3.0 



/^B in the dynamical calculations, respectively, being in 
good agreement with the experimental value 3.2 /^b-^^ 
Calculated Curie temperature is 2020 K (2070 K) in 
the 2nd-order dynamical calculations (the static approx- 
imation). They are much smaller than the Hartree-Fock 
value 12200 K, but still twice as large as the experimen- 
tal value (1040 K).'*" The present results are compara- 
ble to the QMC result of calculations without transverse 
spin fluctuations (1900 K).^^ The reduction of Tq due 
to dynamical corrections is 50 K, which is rather small. 
Dynamical effects in general reduce the magnetic energy, 
but also reduce the magnetic entropy of the static ap- 
proximation. Both effects are competitive to each other, 
resulting in the reduction of Tq by 50 K in the case of 
Fe. 

The magnetization increases with decreasing tempera- 
ture, and reach the Hartree-Fock value 2.61 y^e at T = 
K in the static approximation. The latter is overesti- 
mated as compared with the experimental value 2.216 
/iB-^"'^ The 2nd-order dynamical CPA calculations yield 
M — 2.59 /iB (extrapolated value); the calculations 
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Fig. 7. Calculated magnetization (M), inverse susceptibility 
(x^^)i ^iid amplitude of local moment ((m^)^ ") as a function of 
temperature (T) in Fe. The dynamical results are shown by the 
solid curves. Results in the static approximation are shown by 
dotted curves. Magnetization calculated by the DMFT without 
transverse spin fluctuations^® is also shown by open squares. Ex- 
perimental data of magnetization''^ are shown by +. Note that 
the absolute values of the DMFT magnetization are not given 
in Ref. 29. Thus they are plotted here by assuming that the 
extrapolated value at T = agrees with the experimental one. 



hardly reduce the ground-state magnetization as seen in 
Fig. 7. One has to take into account the higher-order 
electron-electron scattering effects as found in the low- 
density approximation^^ to reduce the magnetization. 
The amplitude of local magnetic moment was calculated 
by means of eq. (84). The results are plotted in the 
same figure. Because of the strong Coulomb interaction, 
it hardly changes with increasing temperature. The dy- 
namical fluctuations enhance the amplitude y {"fn ) by 
1%, and reduce the d charge fluctuations \J {{5hdY) by 
5% at T/Tc = 1.2. 
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Fig. 8. Calculated DOS in the paramagnetic Ni. Solid curve: 2nd- 
order dynamical CPA, dashed curve: static approximation, dot- 
ted curve: XPS data.* 
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Fig. 9. Calculated DOS in the ferromagnetic Ni. Solid curves : 
spin-polarized d DOS in the 2nd-order dynamical CPA, dashed 
curves : spin-polarized d DOS in the static approximation, and 
dotted curve: total DOS in the 2nd-order dynamical CPA. 



We have also calculated the electronic and magnetic 
properties of the fee Ni at finite temperatures. Figure 
8 shows the DOS in the paramagnetic state. In the 
static approximation, the details of the structure are 
smeared by thermal spin fluctuations and the d band 
width is broadened by about 0.1 Ry. The dynamical ef- 
fects suppress the thermal spin fluctuations and develop 
the quasiparticle states. Reduction of the quasiparticle 
band width is 17% as compared with that of the static 
approximation. Furthermore we find the satellite peak at 
Lu =- —0.45 Ry. These results explain well the XPS data® 
as shown in Fig. 8. 

Below Tc, the peak of the down-spin band is on the 
Fermi level, as shown in Fig. 9. On the other hand, the 
top of the up-spin d band is away from the Fermi level, 
so that the peak is weakened due to larger damping of 
the quasiparticle states. The satellite peak for the down- 
spin band disappears due to strong incoherent scatterings 
around lo = —0.35 Ry, while the satellite peak for the up- 
spin band is enhanced sX lo — —0.45 Ry. 

The effective potential for Ni shows a single minimum 
structure as shown in Fig. 10. The minimum position 
shifts to the origin with increasing temperature. This 
should be contrasted to the case of Fe, in which the ef- 
fective potential has a double minimum structure even 
above Tq as shown in Fig. 5, and the paramagnetic state 
is realized by changing the energy difference between the 
two minima. The dynamical potential -Edyn(C) in Ni has 
a butterfly structure as in the case of Fe, but it is highly 
asymmetric along the z axis in the ferromagnetic state so 
that considerable reduction of the magnetization due to 
dynamical corrections occurs. We flnd the reduction of 
\/{il) by 5.0%, and the enhancement of \/ {S,\) by 1.5% 
at T/Tc = 1.3. 

The magnetic moment and the inverse susceptibility 
calculated from the effective potential are presented in 
Fig. 11 as a function of temperature. The susceptibil- 
ity follows the Curie- Weiss law. Both the static and dy- 
namical calculations yield the effective Bohr magneton 
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Fig. 10. Effective potential in the ferromagnetic Ni at T/Tq 
0.47 on the ^a^-^z plane. 
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Fig. 11. Magnetization, inverse susceptibility, and amplitude of 
local moment as a function of temperature in Ni. The dynamical 
results are shown by the solid curves, while the results in the 
static calculation are shown by dotted curves. 



number 1.2 /ie, which should be compared with the ex- 
perimental value 1.6 /ie-*'^ Calculated Curie temperature 
in Ni is 1260 K (1420 K) in the 2nd-order dynamical 
(static) calculations. These values are much smaller than 
the Hartree-Fock value 4950 K, but are still twice as large 
as the experimental value 630 K.^° 

The magnetization increases with decreasing temper- 
ature below Tq. Extrapolated value at T = is 0.67 
(0.71) /iB in the 2nd-order dynamical (static) calcula- 
tions. These values are considerably larger than the ex- 
perimental one (0.62 /ie)-'*^ The amplitude of Ni local 
moment slightly increases with increasing temperature 
and hardly shows anomaly at Tq . The 2nd order dynam- 
ical corrections to the amplitude of local moment and the 
local charge fluctuations are small; the enhancement of 

rh?) is only 1.8 % and the reduction of the d charge 

fluctuations ^JW^df) is 5.9 % at T/Tc ^ 1.3. 



6. Summary 

We have developed the dynamical CPA on the basis 
of the LDA-I-TB-LMTO Hamiltonian towards realistic 
calculations of the itinerant electron system. The the- 
ory is a direct extension of the single-site theory devel- 
oped by Cyrot, Hubbard, Hasegawa, and Kakehashi, to 
the degenerate band case. It is based on the functional 
integral method which transforms the interacting elec- 
tron system into an independent electron system with 
time-dependent random charge and exchange potentials. 
Using the method, we have taken into account the spin 
fluctuations as well as charge fluctuations in the degen- 
erate band system. We have then introduced an effective 
medium YiL„{iojn), and derived the self-consistent dy- 
namical CPA equation for the medium, making use of a 
single-site approximation. 

We adopted the harmonic approximation (HA) to 
treat the functional integrals in the dynamical CPA. The 
HA describes the dynamical effects from the weak- to the 
strong- Coulomb interaction regime. The approximation 
allows us to obtain analytical expressions of the physi- 
cal quantities, and takes into account the dynamical cor- 
rections successively starting from a high-temperature 
approximation {i.e., the static approximation). We can 
calculate the excitation spectra as well as the thermody- 
namic quantities even at low temperatures using the HA 
because we obtained their analytic expressions. 

We have investigated the dynamical effects in Fe and 
Ni within the 2nd-order dynamical CPA, and have shown 
that the 2nd-order dynamical corrections much improve 
the single-particle excitation spectra in these systems. 
The static approximation broadens the DOS due to ther- 
mal spin fluctuations at finite temperatures. The dynam- 
ical effects suppress the thermal spin fluctuations and 
create the quasiparticle states with narrow band width 
near the Fermi level. Furthermore, the correlations cre- 
ate the satellite peak at 6 eV below the Fermi level in 
both Fe and Ni. The XPS data in the paramagnetic Ni 
is well explained by the present theory. 

We verified that the dynamical CPA yields the Curie- 
Weiss susceptibilities. Calculated effective Bohr magne- 
ton numbers, 3.0 /kb for Fe and 1.2 /xb for Ni, explain the 
experimental data quantitatively or semiquantitatively. 
Calculated Curie temperatures, 2020 K for Fe and 1260 K 
for Ni, are however overestimated by a factor of two. Ex- 
trapolated values of the ground state magnetization, 2.59 
/xb for Fe and 0.67 jib for Ni, are also overestimated con- 
siderably as compared with the experimental ones (2.22 
/zb for Fe and 0.62 ^b for Ni). 

We found that the static approximation provides us 
with a good starting point to calculate finite-temperature 
magnetic properties of Fe and Ni, but the dynamical cal- 
culations to go beyond the static approximation have 
been limited to the second-order dynamical CPA in the 
present work. Overestimate of the ground-state magne- 
tization and the Curie temperature should be reduced 
by taking into account the higher-order dynamical fluc- 
tuations. Further improvements of the dynamical CPA 
theory are left for future investigations. 
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We calculate in this appendix the coefficients 
\a \i'^' ^' '^) i^ ^^*^ n-th order expansion of the deter- 
minant D^{k, m) with respect to the dynamical potential 
Va(v,m). 

Let us rewrite D^{k,7n) defined by eq. (57) as follows 
by making use of the Laplace expansion. 
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It should be noted that cq. (A-1) is calcu- 
lated from a set (-D10, -D20, -D21, -D22, -D23, -D24)- 
Thus wc define D^") by *£)(") 

(£'nO, -Dn+10, -Dn+11, -Dn+12, -Dn+13, -Dn+14)- By mak- 
ing use of the Laplace expansion, we can derive a 
recursion relation as follows. 

("'+4"^)£)("+2) . (A.12) 
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In the above equations, we have omitted the suffixes 
v^ k, m for simplicity, and |a'^"''| denotes the determinant 
of the 2x2 matrix ak{v,m). {Dna} at the r.h.s. of eqs. /„(n) 
(A-l-A-5) arc defined by 
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cj^ , and C2 matrices are of the 0-th order, the first or- 
der, and the second order with respect to the dynamical 
potential v, respectively. 

Using the relation Cg = cq and DqoO ~ 1, wc obtain 
the relation. 



Here a'"-* {a^"') stands for a„iy+fc(z^, tti) (a_„i/+fe(i^, m)). 
da , ba , tta , and 5q, are defined by Oq = a'"^ , 6q = 
On = a^^"' bn =1, and for a = 1 ^ 4, 
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Here*£;2 = (1,1,0,0,0,0). 

Substituting cq. (A-12) into eq. (A-15), and using the 
recursion relation successively, we reach the expansion of 
D^ ' with respect to the dynamical potential v. 
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where ii ■ ■ -ik take a value 1 or 2. Here A = (1 + i)/2, and 

In the same way, we obtain the expansion of D^'^ as ^^ ^ ^^« (^ _ ^)^^« (^) (a = 0, a;, y, z) , (A.25) 

at^=g2in-iy)glin)±glin-iy)g2in) ia^x,y,z),{A-26) 
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Here ci"^ and ^"^ are defined^ by ci and C2 in which ^f (x,y,z). The static Green functions gl{n) {a = 

{a(")} have been replaced by {«(")}. 0,a;,y, z) are defined by ~gLaa'{n) (see eq. (5f)) as 

Substituting eqs. (A-16) and (A-17) into eq. (A-f ), we 

obtain the expansion of D with respect to dynamical ~ , ^ o/ \r , V~^ «/ \/ ^ /a oo^ 

potentials 5L^^'(n) == 5/(n)(5^^. + ^^ g^(n)(fj„)^^' . (A-28) 
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Appendix B: Calculation of the Gaussian aver- 

xv„^{v,m)v^^{~iy,m)D\''J^^{iy,k,m). (A- 18) age of dynamical potentials 

Note that «„ and 7„ take 0, x, y, and z. Tir 1 1 ^ 1 iu /-( • r iu ^u 

" '" «(.\ ' ^' We calculate here the Gaussian average of the n-th 

The first few terms of D^^'J_^^ {1^, k, m) are expressed as ^^^^^ products of dynamical potentials, 
follows. 

L>(°)(iy,fc,m) = 1 , (A-19) 

OC 

i)£l) (i^, fc, m) = - ^ y^ ag-X"^-, m, nv + k)„„, (A-20) 
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Here ciajivTm^n) is defined by 

[(1 + Oicr^ + 02cry + 030-2) h{m, n-v, n)] ^^ . (A-22) 
Oi, O2, O3, and /i in eq. (A-22) are 4x4 matrices defined p;+i 

by ~ 



Here integers {n(?Ti)} satisfy the constraint ^^ n{m) — 

n. CMMM stands for Emm' C(^^)^mm'Cm'('^)• 
^'Q,(^/, to.) are given by eqs. (61) and (62). 

The average is calculated from a generating function 
I{s,t) as follows. 
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The latter is obtained as follows, 
r 1 v^ 



I{s, t) = exp 
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Here the Coulomb interactions C^^/ are defined by eq. 
(66). 

By differentiating I{s, t) with respect to Sma 
{t„ij), we have a new factor (2/3)"^ X]„ C'mn^na 
((2/3)"^E„Sn7Cnm)- When we take the 2n-th 
derivative of I{s, t) with respect to the variable 
(smai,tm7i, • • • , Sma„ , im7„ ) , we have a 2n-th order 
polynomial times I{s,t). When we put Sma = and 
tma = in the derivative, we have /(s — 0,t ^ 0) = 1, 
and only the 0-th order terms of the polynomial remain. 
The latters were created by taking a derivative of the 
factor (2/3)-iE„*naQ™ or i2P)-'j:^,CZ^,t^'j with 
respect to the variable conjugate to Sm-yi or tm'ai- A 
created constant (l/2/3)C^^,(5a.-y^. may be indicated by 
a contraction Smaitm'-yj- Then the 0-th order terms, and 
therefore the Gaussian average (B-1) should be given by 
the sum over all possible products of contractions. 
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Here the permutation P is taken with respect to the n 
elements {(fc, m) | fc = 1, • • • , n{m); m = 1,- ■ ■ ,21 + 1}; 
P{{k,m)} ~ {(fcp, TTip)}. Application of the formula 
(B-5) to cq. (63) yields eq. (65) in §4: 
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